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. abstract 

CD ■ 

Ph ' A charge- free, point particle of infinitesimal mass orbiting a Kerr black hole is known to move along a geodesic. When the 

qq , particle has a finite mass or charge, it emits radiation which carries away orbital energy and angular momentum, and the orbit 
deviates from a geodesic. 

In this paper we assume that the deviation is small and show that the half-advanced minus half-retarded field surprisingly 
provides the correct radiation reaction force, in a time-averaged sense, and determines the orbit of the particle. 

> ' 

in 
t> 

O ■ I. INTRODUCTION 

(N 

Binary systems of solar-mass compact objects and supermassive black holes at galactic nuclei are expected to be 
important sources of gravitational waves. In order to detect these gravitational waves, a project to construct a space- 
based detector, LISA, is underway [1]. The detection of such gravitational waves will reveal fundamental information 
\ about gravitational theory, and will provide insights into conditions at the centers of distant galaxies. It may be pos- 
q-i sible to read out the detailed geometrical structure of a supermassive black hole encoded in the incoming gravitational 
waves, with which we can test gravitational theory in the strong gravitational region, and our understanding of black 
K 1 holes. Not only for extracting such information from detected gravitational waves, but also for a more efficient detector 
rr| design of the presently on-going project, it is an urgent problem to calculate as precise as possible the gravitational 
^ , wave signal expected from such a binary system. 

Because of the extreme mass ratio, such a binary system can be treated by a perturbation formalism. We treat the 
supermassive black hole as a background, and treat the solar-mass compact object as a source of metric perturbation. 
It was shown t that, when the spatial volumn of the solar-mass compact object is smaller than the background curvature 
scale, one can approximately use a point source 

p ,, = ( , S w (x-z(T))dz»dz v 



\J—g dr dr 

where z M (r) is the orbit of the object, and r is proper time. 

By the uniqueness theorem of black holes, we can assume that the background black hole is described by the 
Kerr geometry. In this case, given an orbit, there is a formalism to calculate the metric perturbation at infinity [3]. 
Therefore, the remaining problem is how to calculate the orbit z^(t) of the source. 

In the massless limit, the particle moves along a geodesic. When the particle has a mass, it becomes a source of 
gravitational radiation, which carries away orbital energy and angular momentum. Thus the orbit deviates from a 
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T In the framework of metric perturbations, it is not a simple problem to define a point particle. When we take the zero-volumn 
limit, the perturbations become divergent around the particle and the perturbation scheme becomes invalid. In Ref. [2], by 
using a matched asymptotic technique consistent with the perturbation formalism, we showed that the use of a point particle 
is still valid to induce a correct metric pertrubation. 
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geodesic. The deviation from a geodesic can be derived fully by solving an equation of motion with the so-called 
self-force. A general calculation scheme for the gravitational self-force was proposed by Ref. [4]. A similar situation 
happens to a particle with small scalar/electromagnetic charge, and the scalar and electromagnetic self-forces were 
proposed in Ref. [5] . In this paper, starting from a discussion of a symmetry property of the self-force, we propose a 
method to calculate the orbital evolution under scalar/electromagnetic/gravitational radiation reaction. We assume 
that the effect of the radiation reaction is weak, and consider the leading correction to the orbital evolution. 

It is commonly assumed that the orbit evolves in an adiabatic manner, namely, the orbit evolves slowly in its phase 
space. A true geodesic around a Kerr black hole is characterized by the energy E, the angular momentum L and the 
Carter constant . Numerous investigations have been made to calculate the radiation reaction effect on energy E 
and angular momentum L [3] by analyzing the asymptotic gravitational waves at infinity and the horizon. However, 
the calculation of the radiation reaction effect on the Carter constant K is so far an unsolved problem. 

A geodesic equation is a set of second order differential equations of 4 functions {z a }. With a proper time r as 
the parameter characterizing the orbit, we have 7 integral constants, z a (r = 0) and dz a /cLt{t ~ 0), three of which 
are related to E, L and K . In Sec. II, we introduce a specific symmetry property of families of geodesies in the Kerr 
geometry. Using this symmetry, we discuss an important property of the self-force induced by a geodesic, and prove 
that the radiation reaction to the energy E, the angular momentum L and the Carter constant K can be derived by 
use of a radiative Green function (a half-retarded-minus- half-advanced Green function). However, this is not the end 
of the story. We find that the orbit does not evolve in a strictly adiabatic manner in general as the energy E, the 
angular momentum L and the Carter constant K vary on a short time scale. Besides, it is not trivial that the rest of 
constants do not evolve by the self-force. 

In order to complete the calculation of the orbit, we perturbatively integrate the orbital equation with the self-force 
in Sec. III. We figure out the part which secularly evolves by the self-force, and define the 'adiabatic evolution' of 
the orbit in the way we can approximately calculate the orbital evolution. By showing the validity of the adiabatic 
approximation, we propose our new formalism as a conventional tool to predict the gravitational waves detected by 
the future LISA project. 

In this paper, we adopt Boyer-Lindquist coordinates {t, r, 6>, </>}, and M, a are the mass and the spin coefficient of 
the black hole respectively. 



The purpose of this section is to show a simple method to calculate the self-force. The core idea of this new proposal 
is to use a symmetry of the background spacetime together with the whole familiy of geodesies. 

In Subscc.IIA, we first discuss how we can define the family of geodesies. Since we are interested in a particle 
motion as a target of gravitational wave observation, we only consider geodesies rotating around a Kerr black hole 
which neither fall into the horizon nor go to infinity. 

Subsec.II B discusses some symmetry properties of the Kerr spacetime. We apply these symmetry transformations 
to the familiy of geodesies and to the self-force induced on a geodesic in Subsec.II C. Using the result of these 
transformation properties, we discuss a general expression of the evolution equations of the energy, angular momentum, 
and Carter 'constant' in Subsec.II D, and find that a part of the evolution equation can be evaluated by using the 
radiative Green function, which has a great computational advantage. 

In Subsec.II E, we give some comments about practical issues seriously discussed in the self-force problem. 



II. SELF-FORCE AND SYMMETRY 



A. Geodesies around a Kerr black hole 



A general geodesic satisfies the equations 




(2.1) 



(2.2) 



dt 

dX 

# 

dX 



(T, 2 E - 2aMrL) , 

[2aMrE + (p 2 - 2Mr)Lcosec 2 



0] 



(2.3) 



(2.4) 



*We adopt the definition of the Carter constant in Ref. [6]. 
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p 2 = r 2 + a 2 cos 2 9, A = r 2 - 2Mr + a 2 , (2.5) 
S 2 = (r 2 + a 2 )p 2 + 2a 2 Mr sin 2 9. (2.6) 

Here we use A as an orbital parameter related with the proper time r by r = J Q A dAp 2 . 

We consider the case that the radial motion is bounded in n < r < r 2 . One can integrate (2.1), and we have a 
radial periodic solution with respect to A with the period 

A r = 2 / dr ■ . (2.7) 

J ri ^[(r 2 + a 2 )E - aLf - A(r 2 + K) 

We write the periodic solution of the radial equation as 

r(A) = R(E, L,K;X- X r ) = R(E, L, K; A - A r + A r ) , (2.8) 

where we set R(E, L, K; 0) = ri (or, equivalently, R(E, L, K; A r /2) = r 2 ). A r is the integral constant of the first 
differential equation (2.1), and we have a periodicity in A r , namely, A r — > A r + nX r does not change the orbit, where 
n is an arbitrary integer. Because (2.1) is invariant under A — > —A, i?(-E, L, X; A) becomes a symmetric function with 
respect to A as 

R(E,L,K;X) = R(E,L,K;-X). (2.9) 

When the orbit is radially bounded, the 9 motion oscillates symmetrically around 8 = tt/2 in a domain < 9\ < 
9 < 7r — Q\ < 7r [6]. Similarily to the radial motion, we have a 8 periodic solution with respect to A by the period 



= 4 / ' d9 - - = . = t L== - - == . (2.10) 



r Tr/2 , 

^-(aE sin 9- Lcosccfl) 2 - a 2 cos 2 9 + K 
We define the solution of (2.2) as 

9(\) = Q{E, L,K;X- X e ) = Q(E, L,K;X-X e + X e ) , (2.11) 

where we set 6(^, L, K; 0) = 9\ (or, equivalently, <d(E, L, K; X$/2) = n — 9\). Xe is the integral constant in solving 
(2.2), and we have a periodicity in Xg as in the radial equation (2.8). The symmetry property of (2.11) becomes 

e(E,L,K;X) =B(E,L,K;-X). (2.12) 

We write the solutions of (2.3) and (2.4) as 

t(X)=T(E,L,K,X r ,X ;X)+i, (2.13) 

T(E,L,K,X r ,X e ;X)= — (S 2 S - 2aMrL) , (2.14) 

0(X)=<t>(E,L,KXAe;X) + 0, (2.15) 

&(E, L, K, X r , Xg; X) = — [2aMrE + (p 2 - 2Mr)Lcoscc 2 9] , (2.16) 
Jo A 

where i and (j> are the integral constants. We have addition formulae of T(E, L, K, X r , Xg; X) and &(E, L, K, X r , Xg;X) 
as 

T(E, L, K, X r , Xg;X) = T(E, L, K, X r — X X ,X$ — X x ; X - X x ) + T(E, L, K, X r , Xg;X x ) , (2.17) 
&(E, L, K, X r , Xg;X) = *(£7, L, K, X r - X x ,Xg - X x ; X - X x ) + *(£7, L, K, X r , Xg; X x ) . (2.18) 

By using (2.9) and (2.12), the symmetry properties of (2.13) and (2.15) become 

T(E, L, K, Ar, Xg; A) = —T(E, L, K, -X r , -Xg; -A) , <P(E, L, K, X r , Xg; A) = -$(£, L, K, -X r , -Xg; -A) . (2.19) 
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B. t- and ^-translation and Geodesic Preserving Symmetry 



As is well-known, a Kerr spacetime has t- and (^-translation symmetry. Applying the coordinate transformation as 

t = t' + t s ,r = r',6 = 6', <\> = </>' + 4> s , (2.20) 

we have a Kerr metric of the same mass and spin parameter with the new coordinates {F , r' , 8' , </>'}. 

We consider to apply this transformation to a geodesic with a A-translation, A' = A — A^, where X x is an arbitrary 
constant. Using (2.17) and (2.18), we have another geodesic in a new coordinate system as 

r'(X') = R(E, L, K; A' - X' r ) , 6'{X') = Q{E, L, K; A' - X' e ) , (2.21) 
F(X') = T(E, L, K, A;, A^; A') + F , 0'(A') = <&(£, L, K, X' r , X' e ; A') + 4>' , (2.22) 

where the orbital constants, {X' r , X'g, F, </>'}, become 

X' r = A r — X x , X' e = Xg — X x , (2.23) 
F =i+T(E,L,K,X r ,X e ;X x ) -t s , <j>' = (j> + ®(E, L, K, X r , Xg; X x ) — 4> s . (2.24) 

One sees that the orbital constants, E, L and K, are invarinat under this transformation. We can set F, 4>' and either 
X' r or X'g arbitrarily by an appropriate choice of t s , <j) s and X x . 

One cannot set both A^, and X' e zero at the same time since we have only one constant A^ to fix. However, because 
of the periodicity of the radial function (2.8), one can replace A r and Xg by numbers congruent to A r and Xg modulo 
A r and Xg respectively, i.e. A r + n r X r and Xg + ngXg where n r and ng are arbitrary integers. Using this freedom, we 
set X x = X r — n r X r , then we obtain A^, = and X' g = Xg — X r + n r X r + ngXg. When the ratio of A r and Xg is irrational, 
there is a choice of n r and ng with which |A' e | become infinitestimally small, and a geodesic is characterized only by 

E, L and K . In the following, we assume that the ratio of A r and Xg is irrational though we do not set A r and Xg 
zero for the latter convcicncc unless stated. 

Using this transformation property, one can prove a useful formula of a, scalar function geometrically defined along 

a geodesic. As the geodesic is characterzed by 7 constants, {E, L, K, X r , Xg, t, <ft}, we write the scalar function as 
f(E,L,K,X r ,Xg,i,4>o;X). We assume that the function is invariant under t- and (^-translation, then the function is 
indcpend on i and (f>. Since the function is periodic with respect to A r and Xg, one can expand the function with 
discrete Fourier series, e ~ l2m7l K / K-i2-rnr\ e / \e _ gy applying the A-translation with X x = A, we finally have 



f(E, L, K, X r , X e ,i, X) = £ (E, L, K) exp 

We next consider the symmetry as 



X — A r A — A^ 
m — : h n- 



Xg 



(2.25) 



t = -t', r = r',6 = 6',<j> = -<f>' . (2.26) 

By this coordinate transformation, we recover the same line element with the coordinates {t 1 , r', 8' , </>'}. 

We consider the transformation of a geodesic by this symmetry. Since we change the time direction, we transform 
the orbital parameter as A' = —A. Using (2.9), (2.12) and (2.19), a geodesic is transformed to a new geodesic as 

r'(X') = R(E, L, K; X' - A^) , O'(X') = B(E, L, K; X' - X' e ) , (2.27) 
t'(X') = T(E, L, K, X' r , X' g ; A') + F , 4>'{X') = $(E, L, K, X' r , X' g ; A') + $ , (2.28) 

where the orbital constants, {X' r ,X' g ,P,(f>'}, become 

A; = -A r , X'g = -Xg , (2.29) 
F = -t, 4>' = -4>. (2.30) 

Since E, L and K are invariant, by using the t- and (f>- translation symmetry (2.20) in an appropriate manner, the 
new geodesic becomes equal to the original one. For this reason, we call this Geodesic Preserving Symmetry (hereafter 
GPS) transformation. 
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C. Green Function and Self- force 

We denote a various scalar/electromagnetic/gravitational Green function by 

G(x, z) , scalar , 

G(x,z) — { G afJ _(x 1 z)dx a dz^ , electromagnetism , (2.31) 

G a fi ^(x, z)dx a dx /3 dz fJ- dz u , linear gravity . 

and consider t- and ^-translation and GPS transformation property of it. 

Using a symmetric scalar/electromagnetic/gravitational Green function Q sym -(x, z), a retarded and advanced 
scalar/electromagnetic/gravitational Green function, Q ret -/ adv -(x,z) can be written as 

g ret -(x, z) = 26[^{x), z\g s y m - (x, z) , (2.32) 

g adv \x, z) = 20[z, V(x)]G svm - (x, z) , (2.33) 

where £(x) is an arbitrary space-like hypersurface containing x, and 0[£(x), z) = 1 — 0[z, S(x)] is equal to 1 when z 
lies in the past of S(x) and vanishes otherwise. The symmetric Green function is invariant under t- and (^-translation 
and GPS transformation because, in its Hadamard construction [4,5], it is described only by geometrically defined 
bi-tensors invariant under t- and (^-translation and GPS transformation. 

Under t- and (^-translation (2.20), the factor 0[E(x), z] is also invariant and we have 

g ret -(x',z') = g ret -(x, z ), g adv -( x ',z') = g adv -(x, z ). (2.34) 

On the other hand, GPS transformation (2.26) changes the direction of the time and the factor 0[E(x), z] transformes 
as 

0[T,(x'),z'} - 6>[z,S(x)], 6[z',Y,(x')} = 0\S(x),z]. (2.35) 

Thence, by GPS transformation (2.26), a retarded and advanced Green function are transformed to be an advanced 
and retarded Green function respectively as 

g ret -{x',z') = g adv -(x,z), g adv -{x',z') = g ret -(x,z). (2.36) 

We next consider the scalar/elcctromagnetic/gravitational self- force acting on the particle. Because the field induced 
by a point particle diverges along the orbit, we need a regularization calculation to derive the self-force [4,5]. Based 
on the Green function method in calculating the field, an elegant method of regularization was proposed [7], in which 
the self-force can be directly derived from the field calculated by the so-called R-part of a retarded Green function. 
The R-part of a retarded and advanced Green function g R - ret -/ R - adv -( x , z) is schematically defined as 

g R - ret -/ R - adv -{ x ,z) = g ret -/ adv -(x, z )-g s (x,z), (2.37) 

where g s (x,z) is the so-called S-part [7]. It is important to note that the R-part of the half-retarded-minus-half- 
advance Green function becomes a radiative Green function g rad -(x, z) as 

\ (g R - ret - ( x ,z)- g R - adv - (x, z)) = i (g ret - (x, z) - g adv - (x, z )) = g rad - (x, z) . (2.38) 

Similar to the symmetric Green function, the S-part is defined by geometric bi-tensors and is invariant both by t- 
and (/(-translation and GPS transformation. Thence the R-part of the retarded and advanced Green function are still 
invariant under t- and (^-translation, and, by GPS transformation, the R-part becomes 

g R - ret -(x',z r ) = g R - adv -(x, z ), g R - adv -(x',z') = g R - ret -(x, z ). (2.39) 

The scalar/electromagnetic/gravitational self- force is schemetically described as 

F R-ret./R-adv.^ = Um pJ^R-ret./ R-adv.^ x ^ ^.40) 

x — >z(t) 



iR-ret./R 



- adV -( X ) = J dTG R-ret./R-adv. ^ z ( T ^ S ( z ( T fi ; ( 2 .41) 



where <f> R ret -/ R adv - is the R-part of a scalar/electromagnetic/gravitational potential using the R-part of a re- 
tarded/advanced Green function, and we note R ~ ret -/ R - adv - to the self-force to emphasize that it is derived using the 



5 



R-part of the retarded/advanced Green function. S(z(t)) is the source term defined along the orbit. We assume that 
the tensor differential operator i^f^x) is defined to satisfy the normalization condition as F a [](z(r))v a (r) = 0. 

We assume that the self-force is weak and the orbit can be approximated to be a geodesic at each instant of time. 
Using this approximation, we consider to calculate the self-force induced by a geodesic. We write the self-force as a 
vector function of the orbital constants E, L, K, X r , X$, i, (f> and the orbital parameter A as 



F R-ret./R-adv. = pR-ret./ R-adv. ^ ^ R ^ ^ ^ ^ y 

In general, 4-velocity and a self-force transform as 



a , dx' a (dx' a \ (dr\ p 



D (dr\D 



dr' J dr 



dx 13 
dx 7 ^ 



dr 
d? 



(r). 



(2.42) 

(2.43) 
(2.44) 



Using these transformation rules, we apply t- and (^-translation (2.20) and GPS transformation (2.26) to the self-force 
induced by a geodesic (2.42). 

We first consider t- and (^-translation (2.20) with X x = 0. Applying the coordinate transformation, we have 



R— ret. I R— adv. 



(E, L, K, X r , Xe,i+t s ,4> + <j> a ; A) = pR-ret./ R-adv., E ^ L ^ R ^ ~ Xg ^ ^ A) 



(2.45) 



Since the metric is invariant under the transformation, we have F R , ret ^ R adv - = p R ret -/R adv. ^ ^ m s 1 the self-force 



does not depend on t and </>. In the following, we write the self-force as F c 



R— ret./ R— adv. 



(E, L, K, X r , Xg;X). 



Finally we consider GPS transformation of the self-force. Noting (2.39), (2.26) transforms the self-force as 

pR-adv.^ ^ R _^ __~ Xg . _ A) = { _ iyF R-ret. {E ^ ^ R ^ ^. A) ^ 

F*- ret -(E, L, K, -X r , -X e ; -A) = (-l) s F^ adv -(E, L, K, X r , X e ;X), 
where s = 1 for a = t, cfi, and s — for a = r, 9. 



(2.46) 
(2.47) 



D. Evolution of the energy, angular momentum and Carter 'constant' 

As we use A as an orbital parameter, we consider the A derivative of these 'constants' as 



±E 
dX 

—L 

dX 



R—ret./ R—adv. 



R—ret./ R—adv. 



R—ret. I R—adv. 



-p 2 F t R - r ^ R - adv -(E,L,K,X r X;X), 



2 j-iR—ret.J R—adv. 



(E,L, K, X r , Xg; A) . 



2(r 2 + a 2 )p 2 v tpR-ret./ R-adv. + ^r pR-ret./ R-adv.A ^ ^ R ^ ^. y 



(2.48) 
(2.49) 
(2.50) 



Since these are scalar functions defined along a geodesic, we can apply the formula (2.25) to (2.48), (2.49) and (2.50), 
and we obtain 



dX 



R—ret. I R—adv. 



= t R ~ ret - IR ~ advAm - :n) (E, L, K) exp 



. / A - X r X - Xg 

i2tt I m — ; h n 



X r 



Xg 



(2.51) 



where we denote E, L and K by £. By the reality condition, we have (£R-ret./ R-adv. (m,n)y = £R-ret./ R-adv. (-m ,-n) ) 
where * means we take the complex conjugation operation. 
4-velocity of a geodesic transforms by (2.26) as 

v a '{E,L,K,-X r ,-X e ;-X) = (—l) s v a (E, L, K, X r , X$; A) , (2.52) 
where s = for a = t, <f>, and s = 1 for a = r,6. Using (2.47) and (2.52), (2.48), (2.49) and (2.50) transform as 



dX 



R—adv. I R—ret. 



(E, L, K, -A r , -A) = - 



—I 

dX 



R—ret./ R—adv. 



(E, L, K, X r , Xg; A) . 



(2.53) 
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We consider the evolution cuations, (2.48), (2.49) and (2.50), averaged at two orbital points characterized as 
z(E, L, K, A r , Xg, i, </>; A) and z(E, L, K, — A r , — Xg, t 1 , </>'; —A). Using (2.53), one finds the evolution equations are de- 
scribed by the radiative Green function (2.38) instead of the R-part of the retarded/advanced Green function as 



R-ret. 



(E, L, K, X r ,X e ;X) + 



—£ 



R-ret. 



(E, L, K, -A r , -Xg: -A) 



—£ 

dX 



R—adv. 



—£ 



R—adv. 



—e 

d\ 



dX 



~\ rad 



—L 

dX 



-i rad 



-i rad. 



(E, L, K, X r ,X g ;X) + 
(E,L,K, X r ,X g ;X) , 
(E,L,K,X r ,X e ;X) = -p 2 F[ ad (E,L,K, X r ,X g ;X) , 
(E,L,K,X r ,X e ;X) = p 2 F™ d -(E, L, K, X r , Xg; X) , 
(E, L, K, X r , Xg; X) = [2(r 2 + a?)p 2 v t F[ ad - + 2 p 4 v r F™ d 



(E, L, K, -A r , -A ; -A) 



(E,L,K, X r , Xg;X) 



(2.54) 
(2.55) 
(2.56) 
(2.57) 



where is the self-force calculated by the radiative Green function (2.38). 

We write 



dX 



rad. 



(E, L, K, X r , ~X e ;X)=J2 ^ M (E, L, K) exp 



/ A - A r A - X 
i2ir m — h n — 

\ A r Xg 



(2.58) 



then we have 
1 



(2.59) 



— ^£^~ ret - ( m > n ) _|_ (?R—ret. (— m,— n)^ \_ ^gR—adv. (rn.n) _|_ gR—adv. ( — m, — n)^ 

= £rad. (m,n), E ^ K ^ L y 

Thus, half of the expansion coefficients of the evolution equations can be derived by using the radiative Green function. 

We comment that, when the ratio of A r and Xg is irrational, our formula genelarizes the result in Rcf. [8], in which 
it is proven that radiation reaction to the energy and the angular momentum along a whole geodesic can be derived 
by a self-force calculated by a radiative Green function^. By (2.51), the radiation reaction averaged per unit A to the 
energy, angular momentum and Carter 'constant' becomes 

lim [ dX^£ R - ret (E,L,K,X r X;X) = £ R - ret {m (E,L,K) . 

By (2.59), (2.60) agrees with the calculation using the radiative Green function as 

lim J- f\x^£ rad (E,L,K,X r ,Xg;X) = £ R - ret ^(E,L,K) . 

It is notable that the dependence of both A r and Xg vanishes in the end. 



(2.60) 



(2.61) 



6 We also note that our formalism specifies the case that the orbit inducing the self- force can approximated by a geodesic, 
while the formula in Ref. [8] applies to a general orbit in scalar/electromagnetic case. 
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E. Practical issues in calculating the radiative potential 



The primary problem in a Kerr case is we have no conventional method to calculate a electromagnetic potential 
and a metric perturbation induced by a point particle **. The construction of an inhomogeneous solution is unknown 
in general, however, a very simple method to derive homogeneous solutions was proposed [10]. In Ref. [10], it was 
also discussed the derivation of the retarded and advanced Green functions as infinite sums of homogeneous solutions, 
which gives a correct metric perturbation only outside the source. For example, when the particle moves in the radial 
domain r m i n < r < r max , the metric perturbation given in Ref. [10] is correct at r > r max and r < r m i n . 

Though the prescription in Ref. [10] is insufficient for inhomogeneous Green functions, it gives the correct radiative 
Green function since it is just a sum of homogeneous solutions. Suppose we calculate the radiative Green function 
following Ref. [10], it is correct outside the source. However, since it is made as an infinite sum of homogeneous 
solutions by construction, it satisfies the source-free Einstein equations at every radial domain. Thus, it is a correct 
radiative Green function in the whole spacetime. 

III. PERTURB ATIVE EVOLUTION OF AN ORBIT 

To make a definite discussion, we consider that, at A < 0, the particle moves along a geodesic characterized by 
the constants, £ = £q, X r = A r o, X9 = X90, i = to and <fi = <j>Q, and that the self-force begins to act on the orbit 
when A > 0, and deviate from the initial geodesic. In this section, we discuss the deviation of the initial geodesic in 
a perturbative manner. We define /x as the charge or the mass of the orbiting particle normalized by the mass of the 
background black hole, and we consider fx is an infinitcstimally small value as an index of the perturbation. 

In order to see how the orbit evolves by the self-force, we first consider (2.51). We define the deviation of £ from 
the initial value £ by S£. Because we only consider the self- force induced by a geodesic in deriving (2.51), we can 
consistenly derive the evolution of S£ only when 5£ — 0(fi a ), a > 0. The evolution of S£ becomes 

6£(Eo,Lo,K ,XrO,Xeo;X) = X£ R ' ret iE , L , K Q ) + Y,£ H ~ ret - (m ' n \E a , L a , K )e l2 ^ m ^^ +n ^^' , (3.1) 
where the coefficient of the linearly growing term is determined as 

gR-ret. = £R-ret. (0,0) _ (3 2) 

One sees that (3.1) consists of two parts; the secular part and the oscillating part. The oscillating part stays 0(ij}) 
at any A. On the other hand, the secular part grows linearly by A, thus, one can consistently derive S£ only when A 
is 0(/i Q ) , > a > -1. 

Because of this oscillating term, one cannot say that the orbit evolves adiabatically in an exact sense. The oscillating 
part shows the interaction of the orbit and the 'heat bath' of radiation. In the time scale of the order O(/x ), the 
orbit just exchanges the energy, angular momentum with the 'heat bath' and they increase and decrease in the equal 
rate. In the long time scale of the order 0(/i a ) 7 > a > — 1, the energy and angular momentum reserved in the 'heat 
bath' escape into the horizon or away to infinity and the orbital energy and angular momentum tend to flow out to 
the 'heat bath'. Thus, as described by the secular part of (3.1), the orbital energy and angular momentum decrease 
linearly by A. 

Though we do not have an adiabatic evolution in an exact sense, we show that the secular part of the 'constants' £ 
becomes dominant over the oscilating part. If the same thing happens to the rest of 'constants', it seems possible to 
define an 'adiabatic' evolution of the orbit in an approximate sense. For this purpose, we discuss an orbital evolution 
in a perturbative manner. We first consider the evolution of r and 9 coordinates in Subsec.III A, then, t and <f> 
coordinates in Subsec.III B. Subsec.III C gives a plausible definition of an 'adiabatic' evolution of the orbit, which 
approximates the exact orbital evolution by a self-force. Subsec.III D concludes the section with a discussion of a 
gauge dependence of an 'adiabatic' evolution which appears only in gravitational case. 

We define an orbit evolving by a self-force as 

t(X) = t {X) + 6t{X) , r(A) = r (A) + Sr(X) , (3.3) 

9(X) = 6 (X) + S9(X) , 0(A) = 0o(A) + 5<f>(X) , (3.4) 
where {to, ro, 6*o, <fio} is the initial geodesic. For the latter convenience, we define a family of geodesies as 

t(X)=T(E,L,K,X r ,X e ;X)+i, r(A) = R{E,L,K;X — X r ) , (3.5) 

6(X) = ®{E, L,K;X- X r ) , 0(A) = <f>(E, L, K, X r , X e ; A) + 4> , (3.6) 



** Recently, some ideas to calculate the vector potential and the metric perturbation induced by a point source were proposed 
[9]. 
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A. r-motion and ^-motion 

Instead of integrating the equation of motion Dv a /dr = F a , we consider to integrate (2.1) and (2.2) to derive the 
motion of r and 9 coodinates. For a convenience, we write (2.1) and (2.2) as 

dr\ 2 ,„ , fd9^ : 



dx) = v ^ UJ = f/(M) ' (3 ' 7) 

where £ = £ + S£, r = r + Sr and 9 = 9 a + 59. 

Taking the leading order deviation from the initial geodesic, (3.7) become 

2^^ = V(£ ,ro),i5£ i + V(£o,ro),r6r, 2^-^- = U(£ , Bo^SS* + U(£ , r Q ) fi 59 , (3.8) 

where we denote f^5£ l = IeSE + ,f. L 5L + f <K SK. Using 2d 2 r /d\ 2 = V(£ ,r ) tr and 2d 2 6 a /d\ 2 = U(£ o ,0 o ),e, we 
have 

The differential equations (3.9) have singularies because dr /d\ and V(£q, r ) vanish at A = A r + (n/2)A r , and d9 /d\ 
and U{£q,8q) vanish at A = Xg + (n/2)Xg, where n is an integer. 

One must integrate (3.9) such that Sr and 59 are smooth at the singularities. We formally integrate the differential 
equations as 

< 3 '"> 

where we define dVi/dX — Vi/2V and dUi/dX = Uj/2U. Here one must add the integration constants at 

X r + (n + l)A r /2 > A > A r + nX r /2, and ci"' ) at Xg + (n + l)Xg/2 > X > Xg + nXg/2, independently for each integer 
n, such that Sr = and 59 = at A = 0. and 5r and 59 become smooth at the singularities of (3.9). 

In order to determine and together with Vi and Ui, we consider the singular structure of (3.9). We write 

V{£ ,r)=v{£ ,r)(r- ri {£ ))(r 2 {£ )-r) , (3.12) 

u(£ ,e) = u(£ ,e)(e-e 1 (£ ))(n-e 1 (£o)-e) , (3.13) 

where v(£ 0} r) is positive at n < r < r 2 , and u(£o, 9) is positive at 6\ < 9 < tt — 9\. r and 9 n of the initial geodesic 
around the singularities behave as 

n + >i(r 2 - ri)(A - X r - nX r ) 2 + o((X -X r - nX r ) 4 ) , 
ro = { V / ' \ (3-14) 

T2 - \v 2 {r 2 - n)(A - A r - (n + /2)A r ) 2 + 0((A - A r - (n + l/2)A r ) 4 J . 



9 X + iu (7r - 20i)(A - A 9 - nA ) 2 + o((A - A fl - nA e ) 4 ) , 

tt - 0! - i Uo (7r - 2^)(A - A e - (n + /2)A e ) 2 + o((A - A e - (n + l/2)A e ) 4 ) . 



(3.15) 



where ui = i>(£o,?*i), v 2 = v(£o,r 2 ), uq = u(£o,9i) = u(£o,tt — 9\), and n is an integer. Thus, the singular structure 
of \^-{£o,ro) and ^^-(£ ,9 a ) becomes 

XV | -^rj (A - A r - nA r) - + ((A - A r - nA r )°) , 

2T/ °' r ° \^^ 7T (A-A r -(n + /2)A r )- 2 + o((A-A r -(n+l/2)A r )0), 

\-T£^>r ) ^-^-nXg)- 2 + 0{(X-Xg^nXgf), 
2U °' ° \-^^ n (A-A e -(n + /2)A 9 )- 2 + o((A-A 9 -(n + l/2)A e )o). 
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Wc define regularization functions as 



V?(£ ,X-X r o) = i 



2tt 



A r (r 2 - ri) 



r M / 1 1 

Wl I e i27r(A-A r0 )/A r _ ^ e -i27r(A-A r0 )/A r _ ^ 



r 2,' 



1 



y 2 \ gi27r(A-A r0 )/A r _|_ J g-i 2 7r(A-A r0 )/A r _|_ ^ 



(3.18) 



C^(fo,A-A 90 ) = i 



2tt gij 
Ae(r 2 - n) u 



1 



1 



gi 2 7r(A-A()o)/A£) _ J g-i2-nr(A-A eo )/A e _ J 
1 1 



; i2Tr(A-A £ io)/A e _|_ ^ g— i27r(A— A OT ) Ae _|_ ^ 



(3.19) 



One can see that dV^ /dX and dW/dX have the same singular structures as Ui/2U and U^/2U, thus, Vi/2V — dV^ /dX 
and U t i/2U — dU^ /dX are regular and periodic with the period A r and Xg. One can expand these differences with 
discrete Fourier series e l27r ( A -Vo)/A r an( j e i27r(A-A () o)/A £ ^ an( j one can integrate in a regular manner. We formally write 
these integrations as 



Vi - V 1 



(So, A - A r0 ) = (A - X r0 )Vi + ]T VS n \£ )e^" — 

n 

- U\ (£ , X - Xgo) = (X - X eo )Ui + ]T u\ n) (£o)e 



(3.20) 
(3.21) 



where the coefficients of linearly growing terms are 



Vi(£ ) 



X 



dX 



r JO 



2V dX 



(3.22) 



We note that there is an ambiguity in adding integral constants, which the finial result does not depend on. 
Using (3.18), (3.19), (3.20) and (3.21), the first terms of (3.10) and (3.11) can be separeted as 



^ ([K V^SS* + V?6£^ , ^ ([Ui - U}W + U}88^ . 

Now we can see that we successfully regularize the first terms of (3.10) and (3.11) since the divergent behavior of V} 
and Uj and the regular behavior of dr a /dX and d9 /dX cancel each other at the singularities of (3.9) as 



drp 
dX 

d6 { 



dX 



f Jr M + o((A- A r0 -nA r ) 2 ) , 

r 2 ,, + o((A-A r0 -(n + l/2)A r ) 2 ), 

e M + o((A-A 0o -nA 9 ) 2 ) , 
-0 M + 0((A - Xg - (n + 1/2)A 9 ) 2 ) . 



°-U} = 



(3.23) 



(3.24) 



We can further rewrite the first terms of (3.10) and (3.11) using R :i (£o,X — A r o) and Q,i(£ ,X — Ago). By taking 
the £ derivative of the geodesic equation, we have 



d RA IV i d QA 1 U ti 

dX ^ = 27 (£o ' ro) ' dX [m J- 2lT M)) - 



(3.25) 



dX 



dX 



These equations have singularities as (3.9), and must be integrated such that and 0,, are smooth at the singular 
points. Using (3.18), (3.19), (3.20) and (3.21), we have 



dr 



d6 



R,i(£o, X - A r0 ) = ^[[Vi- V?] + V? + c\ v > ) , Q ti (£o, X - A eo ) = ( [Ui - Uj] + Uj + c. 



dX 



rti , rrt , „(«) 



(3.26) 
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where c\ v ^ and are finite integral constants. Using the ambiguity of integral constants in evaluating [Vt — U/] and 
[Ui — U\], we set cf^ = c\ u ^ = 0. and we have the first terms of (3.10) and (3.11) as 

R i (£ ,X-Ko)S£ i , e,i(£ ,A- X eo )S£ l . (3.27) 
Since r- and ^-motion of a geodesic is periodic, we can put 

R{£,X-X r ) = ^fl ( "'(Oe ,2m ^, 6(£, A - A r ) = ^ 6 (n) (£ ) e i2 ™W . (3.28) 

n n 

The ^-derivative of R and becomes 

fi i= V(_ & ^^ R W +fl wl 12 ™^ = -(A-VJ^+Vf^, (3.29) 
V A r X r ' / aA A r ^ 

n x 7 n 

e, 4 = ^-,2™— — e<> + e^j e * = -(a- a*)-— + £e^ e (3.30) 

One can see that R A and 9, are dominated by oscillating parts whose amplitude grows linearly in A. Comparing 
(3.20), (3.21) and (3.26), we find 

Vi(£o) = -^(So), Ui(Eo) = -^(£o). (3.31) 
We next discuss the second terms of (3.10) and (3.11). Using (3.18), (3.19) (3.20) and (3.21), we separete the terms 

as 

-Wl dm - v ' ] i s£ ' + C ixv 'i 5e ')- ,3 - 32) 

d6 ' ' " 



d xy Jo dX ^- u 'W El + h dX ^Tx d£l )- ^ 

The integrands of the first terms in the brackets are regular and periodic, and one can evaluate the integration in 
an usual manner as 

t dX[ V t - Vft-jLsS* = (A ~ 2 Ar ° )2 U e + J2 (( A - Ko)V e (m ' n) + K (m '" } ) e l2 <™^+^) , (3.34) 

[ X dX[U t - Ut\±&(* = (A ~ o Ae0)2 t/ e + V ( (A - Xeo)U^ n) + U^AM mh ^ +nX -^r) . (3.35) 
Jo "A 2 — ' V / 

" u m,n x ' 

The second terms in the brackets could lead to logarithmic divergence at the singular points, r = ri,r 2 and 
9o = #i,7r — 6\. If we have logarithmic divergence, we have no way to have a smooth evolution of the orbit at the 
singularities. Thus, one must constraint on the self-force as 

= r M -^-<5£ 4 (A = X r0 + nX r ) = r 2 , i -^5£ i {X = X r0 + (n + 1/2) K) , (3-36) 

= hl ^-S£ l (X = A eo + nX B ) = e^&S^X = X eo + (n + l/2)X e ) . (3.37) 
d\ d\ 

Using (2.51), we can re-expand as 

hi dX = ^ rl ~ ' 

r2 > i 7\ S£i= E ^r n) ((-l) ro e i2ro7r ^ a -l)e <2 "' r ^ a , (3.39) 

6i t £6£* = V efcrtj^^ie** 1 *^ - 1) , (3.40) 



■'dX^' = 



1\„ —8£ i = ^'"^^"^((-lre^-jT -1), (3.41) 



n 



where £%" n) = nJ l ( m ' n \ £%" n) = {-l) m r 2 ,i£ i{m ' n \ £ { e \' n) = 9 hi £^ m ^ and £^" n) = -(-l) n 9 hi £ i{ - m > n \ Using 
these expansions, we can integrate the second term without logarithmic divergence as 

[ X dX V^S£* = (A-A r0 )yt + VytM*(^+»^) , (3.42) 

JK) m.n 

f X dXU}^-5£ 1 = (A - \ 60 )Ui + V ^("^V 2 ^" 1 ^^"^ 1 ) . (3.43) 
J ° dX 

It is notable that we have only linearly growing terms, thus, (3.34) and (3.35) dominate over (3.42) and (3.43) at 
A - 0(p a ), > a > -1/2. 

Using (3.26), one can rewrite the second terms of (3.10) and (3.11) as 

dr o f X „Ri d r „, d9 f x „8i d 



dX J ifc dX dX J fa d\ 1 ' 

By the constraints on the self-force (3.36) and (3.37), the integrands of (3.44) have no singularity, and we can formally 
write as 

l^j\ 5£i = E(( A -^o)^'" ) +^ m ' rl) )e i2 ^ m ^ a+n ^ 11 ^ , (3.45) 
m^ 5£i = E (( A - ~ X ooM n - n) + &i m n) ) e M™^ iL + n ^ L ) , (3.46) 

dX m.n ^ ' 

where the coefficients of linearly growing terms become 

j^(m.n) ^r.i £R—ret. (m.n) i Q(ni.n) ^6 ,i £R—ret. (m.n) i (3 47) 

A r Xg 

We show how to integrate (3.10) and (3.11) such that 5r and 59 evolve smoothly, and we find ci™' = = 
along this integration procedure. The perturbative evolution of 5r and 59 by the self-force is now interpreted as the 
evolution of the orbital 'constants', 5£ l , 5X r and 5Xg as 

5r = R.M 1 + R - Xr 5X r , (3.48) 
S~ Xr = (A ' 2 Ar ° )2 A r + - X ro)\l m - n) + Xi^Ae a < m ^ +n ^) , (3.49) 

m,n ^ ' 

59 = Q. ;l 5£ l + Q- Xe 5X e , (3.50) 
S X 9 = i^plxe + J2 ((A - A eo )A<r i) + A^ n) ) M mX -^ +nX -^) . (3.51) 

m.n ^ ' 

From (3.47), we have 

X r = = -blA£R-ret. (0,0) ^ = _ ^R-ret. (0,0) i ^ q.V2) 

X r Xe 

Um.n) _ 1 X r.i fR-ret. (m.n) i \(m,n) _ \ ^8 ,i £R- r et. (m.n) i (3 53) 

A (0,0) = ^(0,0) = ^(0,0) + yl ) X (m = @(0,0) = £(0,0) + jjt ^ (3 54) 

X ^ ,n) = 2^f^ + ^) ^ m '" ) = + V * {m ' n) ' = 2n(^ + ^ 9 ° m ' n) = f/e(m ' n) + ^^'^ ' (3 - 55) 

where the second line and the fourth line are evaluated for (m, n) ^ (0, 0), and a[ 0,0 \ A^ 0,0 ' are determined such that 
5X = 5Xg = at A = 0. 

We can see that the linear perturbation holds for A <~ (3(/i Q ), > a > —1/2. The orbital evolution in this regime 
has two parts; the one contributed by 5£, the_other by 5X r and 5Xg. Since both parts grow quadratically by A, it is 
necessary to consider the evolution of A r and Xg for a correct orbital prediction. 
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B. t-motion and 0-motion 

We write (2.3) and (2.4) as 



^■=X(£,r,6), ^=Y{E,r,0). (3.56) 
The leading deviation from the initial geodesic satisfies 

^ - Xi^o, r 0) 6 )5£ 1 + X tr (£o, r ,9 )Sr + X i6 (£ Q , r , 9 )S6 

= (I; + X rJ R ;i + X e e ti )S£ i + X, r R - K d\ r + XjQ^SXe , (3.57) 

^ = r,(£ , ro, 0o)<^ + r r (f , r , #o)<5r + r , 0o)<50 

dA 

= (F, + r r i?, 4 + Y fi Q,i)5£ l + Y r R x SX r + Y Q- X$ 5\e , (3.58) 
where we use (3.48) and (3.50). Contrary to (3.9), (3.57) and (3.58) are regular, and we can integrate by parts as 

St = X£& + X r SX r + XgSXg - ( d\(x~S£ l + X r ^-5X r + Xe^rS'Xg) , (3.59) 

J Q \ dX dX dX J 



— Xi = X t i + X ir R t i + X^Q^i , ~j^X r = X^ r R - Xr , -j^ x 8 = XjQ Xe > (3.60) 



= Y5£ l + Y r SX r + Y SX e - jf dx[f~5£ l + Y r -^5X r + Ye^SXg^j , (3.61) 



= Y l+ Y r R ti + Y e @,i , ^Y r = Y r R x , ^Y e = Y. e Q- Xg . (3.62) 



We note that T :i (£ ,X r o, Xgo; A) and $^(£0, X r o, Xgo; X) satisfy 

^-Ti = X ; + X >r R ti + X fl e,, , - X,i + X >r R ti + X t oQ,i , (3.63) 

dA dA 

thus, from (3.60) and (3.62), we have 

Xi = Ti(£ ,X r0 ,X e o) + c^ , Yi = *,i(£o,Xro,Xeo) + cM, (3.64) 

where and are integral constants. Because the orbital evolution does not depend on and in the end, 
we set them zero in this subsection. In the same manner, we have 

X r = T X , Xo = T- Xo , Y r = , Y B = <S>- Xo . (3.65) 

We define t- and 0-motion of a geodesic as 

T(£,A r ,A e ;A) = AT(£) + ^ (£)e J2 "( m ^ + ™^> ) , (3.66) 

m,n 

X r , X g ;X) = X<Z>{£) + ^ & m > n )(£ y 2 *( m ^m+ n ^m) , (3.67) 

where the linearly growing terms appear since we integrate X(£,R, Q) and Y(£,R, G) which can be expanded by 
discrete Fourier-series e i^m(X-X r )/X r +i2nn(X-X e )/X e _ Us i n g (3 64) and (3.65), we have 

Xi = J2( xx t' U) + X J (m '" ) )e j2 "( m ^ a+ "^F a ) , % = J2( XY l m ' n) + y/ m '" ) )e t2 "( m ^ a+ "^ M ) , (3.68) 

m.n m.n 

X r = J2x^e i2 <^^+^) ?r = j^ Yr(m . n)e ^( m ^ +n ^) (3 69) 

m.n m.n 

X e = J2x^e i2 ^ m ^ +n ^r) , Ye = ]T Y M^( m ^ +n ^) (3 ?0) 
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where the coefficients become 

X| m ' n) = 5 m=n =ot, i - i2w (m^ + n^j , x| m ' n) = T ( /"' n) + i2n (m^f* + n^^j , (3.71) 

Y} m ' n) = S m = n = <f>,i i2* (m^- + n^j , Y^ n) = ^ + z2k (m^ + n^^j , (3.72) 

A r Xg X r Xg 

We consider to evaluate the integration of (3.59) and (3.61). We first note 

*6K = %i4t8S i , 4*^ = ^4^*. ( 3 - 74 ) 
dX # dX dX # dX J 

and the expansion of RHSs are given in (3.42), (3.43) and (3.47). Using these results, the integrands of (3.59) and 
(3.61) become 

^(XX { e m - n) + X^ n) )e 2 ^ mh ^ +nhZ ^> , ^(XY e { " hn) +yJ m <")) e 42T ( m± ^ +n ^F 11 ) , (3.75) 

m.n m,n 

where the coefficients of the linearly growing and oscillating terms become 

j^(m.n) rp gR—ret. (m,n)i y(m,n) ^ gR—ret.(m,n)i ^ JQ^ 

In summary, we have t- and 0-motion as 

5t = T.58 1 + T- K SX r + T- Xe SX g + 5t, (3.77) 
6i = yi'o + Y.^^ + 4 m '" ) )e l27r ( m ^ a+ "^ M ) , (3.78) 

m,n 

Scb = * 4 <5r + 1> x 5X r + <P rXe 5X e + H , (3.79) 

<rl ^ 1 i \ i( m >") i A m - n )\ »27r(m A ~/ r0 +ra A /[ eo ) foon\ 

H= ^^o+l^iMa + <Ph ) e V r 6 . ( 3 - 80 ) 

m.n 

where the part of coefficients of <5i and <5</!> are 

to = -Ti£ R - ret - (°>°) 1 , = 4 £ i? ~ ret - (0 < 0) 4 , (3.81) 

l(m,n) _ * ™ cR-ret.(m,n)i l(m,n) _ 1 i cR-ret. (m.n) i i \ / / n n \ S9\ 

~ 2tt(^ + ' ' ™ ~ 2tt(^- + ^-) ' ' ( . m ' n ^^ ' ^- ^ Al > 

V / V A r \q ) 

As in the case of r- and ^-motion, the linear perturbation holds for A <~ 0(/i Q ), > a > —1/2'', and 58, 5X r , 5Xg, 
5i and 5(f) contribute to the orbital evolution equally. 

C. Adiabatic Evolution 

In Subscc.III A and Subscc.IIIB, we find the evolution of the orbit by the self-force in a perturbative manner. For 
the definite discussion, we assume that the self-force begins to act on the orbit at A = 0, and find that the perturbative 
evolution is correct at A = 0(p, a ), > a > — 1/2. Using this result, we first consider to rewrite the orbital equation 
in a numerically convenient form. 

When the orbit is a geodesic, r- and ^-motions are periodic. Then it is convenient to describe a geodesic with phase 
functions Xr and \e as 



^This time scale is called the dephasing time. [12] 
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R(£, Xr ) =Y,R {n \£V 2mXr ■ ^Xr = , (3.83) 

e(£, X e) = J2e {n) (£)e i2 ™ x °, ^ Xe = xj£)- (3 ' 84) 

We consider the similar type of description when we consider the self-force. One can rewrite (3.48) and (3.50) as 

5r = R t 6£ l + R x Sxr , Ri = Y1 R { f^ nXTQ , R x = Yl *27mi?("V 2 ™ x "° , (3.85) 

n n 

59 = QM l + Q x 6xe , 6, = J2 O^e' 2 ™*" , 6 X = J2 i^ne^e* 2 ™*™ , (3.86) 

n n 

where Xro an d Xeo describe the phase functions of the initial geodesic, and Sx r and 5xe are the deviation by the 
self-force satisfying 

A — A r o A r i , 5X r0 A — Ago \e,i :c i SXeo a on\ 

5xr = r r— St — , 5xe = r t—S£ — . (3.87) 

A r A r A r Ag Ag Ag 

We find that, with the effect of the self-force, the phase functions satisfy 

^ (Xr + S Xr0 ) = A - I _ &£i + F Xr (£ , XrO, XOo) , F Xr = _ J_ ^ Jj(-.n) e «-(-XH,+nx.o) , (3 . 88) 
r r,z r m,n 

^(Xo + Sxeo)= x + [ _ S£i +F Xa (£ , X ro,Xeo), F Xo = -i^eMe 8 "^'). (3.89) 

' l ' m,n 

One can see that (3.87) becomes that of a geodesic when we switch off the self-force, thus, by the renormalized 
perturbation method, we formally obtain the orbital evolution eqution with the self-force as 



= F £ (£, Xr ,Xe) , ^jXr = F Xr (£, X r,Xe) , ^jXe = + F xe (£ , X r, Xe) , (3.90) 



r = V flW (£)e l2 ™ x - , 9 = V e (n) (£)e t2miX( > , = X(£, r, 6) , -^-0 = F (£, r, 0) , (3.91) 

n n 

where F £ is given by (2.48), (2.49) and (2.50). 

The evaluation of (3.90) and (3.91) may be possible by a future investigation of a regularization calculation [11]. 
But, the calculation may be complex and costly unless we have a new breakthrough in the regularization calculation 
strategy. We, therefore, introduce an approximate calculation using only (2.61), of which we already have a well- 
established calculation technique and a number of results [3] . 

From the perturbative result in Subsec.III A and Subsec.IIIB, the orbital evolution is dominantely described as 

5£ = A£*-ret.(o,o) (£q) ; ^ r = - ^ ^£ R ~ ret - (°'°> \£ ) , 5\ e = _^^i£*-™*-(°.°)*(£ ), (3.92) 

2 A r 2 Ag 

5t = ~^Tj R - ret - (0 < 0)t (£ ), S4> = -^J R - ret - {0 ' 0)l (£ ). (3.93) 

We call the calculation using these dominant parts by an approximate adiabatic calculation, and the evolution equation 
becomes 

^ cadi. pR— ret. ( cadi. \ ^ adi. ^ ^ adi. \ /o (\A\ 

dX ( h d\ Xr ~A r (£«*-)' d\ Xt> ~ X e (£^)' ( > 



„adi. 



y2R( n \£ adi -)e i2 ™ x ° di ' , e adt = y2e( n) (£ adi -)e i2nnx e di ' , (3.95) 



— t adu = x(£ adi -,r adi -,e adl -) , — (j) adi - = Y(£ adi -,r adL ,e adi -) . (3.96) 
dX dX 

For a practical use of this approximate adiabatic calculation, we make a rough estimate on how correctly the orbit 
can be predicted by this method. We suppose to divide the whole domain of A into an infinite number of domains 
Xk < X < Xk+i whose interval is 0(/z"), a — > —1/2. Then we may apply our perturbative analysis of the orbital 
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evolution at each domain. We define the orbital 'constants' at A = A^ as {£^ k \ \i k \ Ag fc \ t^ k \ cf)^}. Again, for a 
definite discussion, we assume that the self- force begins to act at A = A . 

By the result of Subscc.III A and Subsec.IIIB, a dominant contribution makes the evolution of these 'constants' as 

£ {k+i) _ £ {k) „ {v 1+a ) , A< fc+1 ) - A« ~ 0(^+ 2a ) , If +1) - \ { e k) ~ 0( M 1+2Q ) , 
fk+l) _ fk) _ 0(iU i+2a) ; ^(fc+1) _ ^(fe) „ (fi 1+2a ) . (3.97) 

After passing by N ~ 0(p ) finite domains, we have 

~ £V+ Q+ ^) , AW ~ fV +2a+/3 ) , A<"> ~ 0(/x 1+2a+/J ) , 
*W ~ 0(/x 1+2a+/3 ) , ~ 0(/i 1+2Q +^) . (3.98) 

On the other hand, the part we ignore for the approximate adiabatic calculation contaminates the evolution as 

5 (£(fc+i) _ f(fc)) „ ( M ) , ,j(A(fc+i) _ AW) ~ 0(/x 1+a ) , <5(A< fc+1) - A< fc) ) ~ 0(/i 1+a ) , 

^ f«) ~ 0(/i 1+a ) , <5(<P +1 ) - <p>) - 0( M 1+tt ) . (3.99) 

We assume that the ignored part affects the evolution like an randam Gaussian noise and estmate the error as 

5tW ~ 0( M 1+a+/3/2 ) , ^ (JV) ~ 0(fi 1+a +^ 2 ) . (3.100) 

Using this estimation, we discuss the predictability of the approximate adiabatic calculate. The error of the rotation 
per time <j)/t is estimated as 



0(p- a - p/2 ) • (3.101) 



For an accurate prediction, we require 5{4>/t) < 1 and we have < —2a — > 1, which is always satisfied since N is just 
an integer. 

The error of the rotation phase <f) is estimated as 

<5</>~ 0(fi 1+a+ P/ 2 ) . (3.102) 

For a correct prediction of the phase, we have > —2a — 2 — > —1. This shows that we have a prediction of the 
rotation phase only at A < 0(/i~ 3 / 2 ). 



D. Gauge Issue 



The gravitational self-force problem has an exceptional difficulty because of the so-called gauge problem [13]. 
The regularization formulation was originally formulated in the harmonic gauge condition [4,7], and we have the 
divergent S-part only in the harmonic gauge. It was pointed out that, when we subtract the S-part from the full 
metric perturbation in the radiation gauge [10], we have a divergent residue [13] because of the divergent gauge 
transformation at the particle location. For this reason, the calculation of the full metric perturbation in the harmonic 
gauge condition, or the calculation of the S-part of the metric perturbation in the radiation gauge become important 
issues in calculating the gravitational self-force in a Kerr background. 

On the other hand, in the approximate adiabatic calculation, all we need is to evaluate (2.61), which is proven to 
agree with the radiation reaction calculation [8] in part. A number of previous works [3] prove there is no divergence 
in the approximate adiabatic calculation. The investigation in Subsec.IID shows that the self-force derived by the 
radiative Green function corresponds to the two-point averaged self-force derived by the R-part of the retarded Green 
function. We consider that, by taking a two point average of the self-force, the divergent S-part vanishes as in (2.38). 
We consider that, in the approximate adiabatic calculation, the same cancellation mechanism happens including 
divergence by the gauge transformation. But there still be an abiguity of a finite gauge choice. Here we prove the 
result by the approximate adiabatic calculation is actually gauge-invariant by showing the gauge invariance of (2.61). 

Using a Killing vector r\ a and a Killing- Yano tensor r\ a p of a Kerr spacetime, conserved quantities along a geodesic 
are described as 

E = Va v a , Q = ri a pv a vP , (3.103) 
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where v a (r) — dz a /dr is 4- velocity and r is proper time of the geodesic. (2.61) can be rewritten as 



lim — — 

T->oo 2T 



lim 



1 



, d r, 1 
dr——E = lim — — 

dr—Q = lim 



dr^Wr))^). 



dr2r; Q/3 (z(T))^(T)^(T) 



T^oo2TJ_ T dr T^oo2TJ_ T 

where F 13 is the self-force vector. 

By a gauge transformation x a — > i a + £ Q , the self- force is transformed as 

F a (r)^ F a (r)+6 s F a (T), 

S s F a (r) = -(^(r) V T(r)f ;/37 ^(r))+i2 Q ^ tf («(r))^(r)r(«(r))t;*(r) 

The increment of (3.104) and (3.105) by this extra term 8^F a (r) becomes 

1 T 



lim 

T^oo 


2T„ 


/ drS 




= lim 

T^oo 


1 

2T 


lim 

T^oo 


2T„ 






= lim 

T^oo 


1 

2T 



(3.104) 
(3.105) 



(3.106) 
(3.107) 



(3.108) 
(3.109) 



where we use Killing equations ?y Q ;/3 7 = r]sR s 7 0a and rj a fj. n s — -q^^s-ya + VaeR e SiP- Since the gauge dependence is 
totally integrated out, the gauge dependence of the approximate adiabatic calculation vanishes by taking T — > oo. 
We comment that the gauge dependence of (3.90) and (3.91) is highly non-trivial and we left it as a future problem. 



IV. SUMMARY 



In this paper, we discuss a method to calculate an orbital evolution by a scalar/clectromagnctic/gravitational self- 
force. We assume that the self-force is weak and that the orbit can be approximated by a geodesic at each instant of 
time, with which one can derive the self-force. We note that the geodesic equation is a set of second-order differential 
equations of four components, and we have 7 integral constants, E, L, K, X , Ai, to an d 4>o- Instead of calculating 
the orbit itself, we derive the equations of these 'constants' by the self-force under this assumption. 

We first consider the evolution of E, L and K since the evolution equations of these 'constants' are directly derived 
by the self-force vector. We exploit the symmetry of a Kerr spacetime together with a family of geodesies, which 
induce the self-force. By applying the symmetry transformation to the self-force vector, we find that the time-averaged 
evolution of E, L and K can be derived by using a radiative Green function, which has a number of technical advantage 
in practice. However, we also find that the orbit does not evolve adiabatically in an exact sense. 

In order to understand the orbital evolution by the self-force, we next consider the orbital equation in a perturbative 
manner. We integrate the orbital equation by a time scale, sufficiently long, but, less than the dephasing time when 
the linear perturbation of the orbit becomes invalid. Since the orbit does not evolve in an adiabatic manner, the 
orbital constants of a geodesic {£ , A r , Xe,i, 4>o} are oscillating by the self-force. However, we could find out secularly 
growing parts which will dominate the orbital evolution. By taking these growing parts only, we define an approximate 
orbital eequation, which we call an approximate adiabatic calculation. We consider that the approximate adiabatic 
calculation is enough implcmcntable by a well-established method since it only uses the radiative Green function. 

We also discuss how approximate an orbital evolution can be obtained by this calculation method. We find that, 
during the time 0{^ z / 2 ), it gives an accurate rotation phase. For example, when 10 solar-mass black hole is 
inspiralling into 10 7 solar-mass supermassive black hole, this corresponds to around 10 9 -rotation period. Though the 
accuracy in predicting a wave form is not clear in our estimate, if only the accurate prediction of the rotation phase 
is important for the future LISA observation, the approximation proposed here may give a sufficient imformation in 
this case. 

Finally we prove that the approximate adiabatic calculation gives a gauge invariant prediction, thus, the result is 
consistent with that by another possible method within the approximation scheme. 
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APPENDIX: CIRCULAR OR EQUATORIAL ORBIT 

When the orbit is circular, or equatorial, we have dro/dX = V(£o,r ) = V(£o,ro), r = 0, or d9 /d\ = U(£o,9 n ) = 
U(£o, Oq)_0 = 0. Then (3.8) becomes trivial and we cannot evaluate the perturbed orbital equation. Here we consider 
the orbital evolution in these special cases by taking the circular, or equatorial limit of a general orbit. Before the 
discussion, we note that, by repeating the symmetry argument in Sec. II, the self-force acting on E, L and K becomes 



dX 



- R—ret./R—adv./rad. 



J2 n £R-ret./R-adv./rad.(0,n)( E ^ LjK y2™-^ ^ c i rcu l a r Orbit , 

J2 m £R-ret./R-adv./rad.(m,0)^ EjL ^ K y2nm^ ^ equatorial Orbit , ( Al ) 

^ £R-ret./R-adv./rad.(0,0)( E ^ ^ K ^ ) circular and equatorial orbit . 



We consider r-motion when we consider a self-force acting on a circular orbit. The perturbation equation (3.8) 
becomes trivial since we have V.i = V. r = along the orbit by using (3.12). Instead of calculating the orbital equation, 
we consider to take the circular limit of (3.48). The first term of (3.48) is given in (3.27) and behaves in a regular 
manner when we take a circular limit. On the other hand, the second term of (3.48) is given in (3.44). By the 
regularization calculation starting from (3.32), the integration is still finite in the circular limit, and dr^/dX vanishes. 
Thus, we have 

Sr = R.M 1 , (A2) 

which means the circular orbit stays circular under the self-force and the orbit is solely determined by E, L and K ^ 
As a result, the orbital equation becomes 

±£ = F £ (£, Xe ), ± X e = -J^ + F XB (£, X r,Xe) , (A3) 

r = R^(£), 6 = J2Q {n) (£)e i2 ™ Xe , ^t = X(£,r,6), ±<f> = Y{£,r,0) , (A4) 

n 

and, under the approximate adiabatic calculation, we have 

^ cadi. cR— ret. I cadi. \ j ,,adi. \ ( A r\ 

dX b ~ b {t > ' dX Xe ~ Xe{£ ad% -) ' ( 5) 

r adi. _ ^(0) (gadi.^ Qadi. _ q(u) ^gadi. ^ e i2imx a e dl ' (A6) 

n 

^ ^adi. j£^gadi. ^adi. Qadi.^ ^ ^adi. -y ^gadi. ^adi. Qadi.^ (A7) 

dX dX 

Similarly, 0-motion can be drived such that the equatorial motion stays equatorial. One can also prove by using 
the symmetry 9 — > —9. Using this symmetry property, when the orbit is equatorial, the self-force satisfies Fg = 0, 
and we have 9 = tt/2, d9/dX = 0. The orbital equation becomes 



dx £-F £ (£, X r), d)y Xr - Xfj{£) 



£ = F £ (£, Xr ), —Xr = ^^ + F Xe (£, X e), (A8) 



r = Y J R (n) {£Y 2lTnXr , 0=\, = X(£,r,n/2), ±<f> = Y(£, r, tt/2) , (A9) 

n 

and, under the approximate adiabatic calculation, we have 

cadi. cR— ret. ( cadi. \ J^_.,adi. ' / a i p»\ 

dX b ~ b {t > ' dX Xr ~ A r (£«<"■) ' (A1U) 

r adi. = ^2 R (n)^adi.y2nn X a r di - ) 6» adl — 7r/2, (All) 
n 

^-t adl - = X(£ adi -,r adi -,Tr/2), ^-4> adl - =Y(£ adi -,r adL ,Ti/2). (A12) 
dX dX 



**The same proof was given in Ref. [14]. 
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When the orbit is circular and equatorial, the oribtal equation becomes 




Fe{£) 



(A13) 



r 



(A14) 



and the orbit evolves adiabatically in an exact sense. 
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